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ABSTRACT 

We show that a purely kinetic approach to the excitation of waves by cosmic rays in 
the vicinity of a shock front leads to predict the appearance of a non-alfvenic fastly 
growing m ode which has the same dispersion relation as that previously found by 



Belli (120041) by treating the plasma in the MHD approximation. The kinetic approach 
allows us to investigate the dependence of the dispersion relation of these waves on 
the microphysics of the current which compensates the cosmic ray flow. We also 
show that a resonant and a non-resonant mode may appear at the same time and one 
of the two may become dominant on the other depending on the conditions in the 
acceleration region. We discuss the role of the unstable modes for magnetic field 
ampUfication and particle acceleration in supernova remnants at different stages of 
the remnant evolution. 
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1 INTRODUCTION 

The problem of magnetic field amplification at shocks is central to the investigation of cosmic ray 
acceleration in supernova remnants. The level of scattering provided by the interstellar medium 
turbulent magnetic field is insufficient to account for cosmic rays with energy above a few GeV, 
so that magnetic field amplification and large scattering rates are required if energies around 
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the knee are to be reached. The chief mechanism which may be responsible for such fields is 
the excitation of streaming i nstabilities (SI) by the sarn e particles which are being accelerated 



Skillingl(ll975h : 



BelK 



1978h : 



Lagage & CesarskyI (|l983al ]bl)). The effect of magnetic field amplifi- 



cation on the maxirnum ene rgy reachable at supernova remnant (SNR) shocks was investigated by 
Lagage & CesarskyI (|l983alJb|) . who reached the conclusion that cosmic rays could be accelerated 
up to energies of order ~ 10"* - 10^ GeV at the beginning of the Sedov phase. This conclusion was 
primarily based on the assumption of Bohm diffusion and a saturation level for the induced turbu- 
lent field 6BIB ~ 1. On the other hand, recent observat ions of the X-ray surface brightness of the 



VolketaL 



(|2005h for a review of results), 



rims of SNRs have shown that dB/B - 100- 1000 (see 
thereby renewing the interest in the mechanism of magnetic field amplification and in establishing 
its saturation level. It is however worth to recall that the interpretation of the X-ray observations 
is not yet unique: the narrow rims observed i n the X-ray sync hrotron emission could be due to 



the damping of the downstream magnetic field 



PohletaL 



(|2005h rather than to severe synchrotron 



loss es of very high energ y electrons, although this interpretation has some serious shortcomings 
(see[ 



Morlino et al 



(|2008f ) for a discussion). 



In thi s con t ext of excitement, due to the implications of these discoveries for the origin of cos- 



mic rays. 



Belli (|2004|) discussed the excitation of modes in a plasma treated in the MHD approx- 



imation and found that a new, purely growing, non-alfvenic mode appears for high acceleration 
efficiencies. The author predicted a saturation of this SI at the level SB/B ~ MAirjvJcY^^ where 77 
is the cosmic ray pressure in units of the kinetic pressure pvj, Vs is the shock speed and Ma = vJva 
is the Alfvenic Mach number. For comparison, standard SI for resonant wave-particle interactions 
leads to expect 6B/B ~ M^^t/^^^. For efficient acceleration rj ~ 1, and typically, for shocks in the 
interstellar medium. Ma ~ 10"^. Therefore Bell's mode leads to 6B/B ~ 300 - 1000 w hile the stan- 
dard SI gives 5B/B ~ 30. It is also useful to notice that the saturation level predicted by 



BelK 



2004) 



is basically independent of the value of the background field, since 6B^/Sn ~ {1/2){vJc)Pcr, 
where Pcr is the cosmic ray pressure at the shock surface. 

The resonant and non-resonant mode have different properties in other respects as well. A key 
feature consists in the different wavelengths that are excited. The resonant mode with the maxi- 
mum growth rate has wavenumber k such that kr^ o = 1, where r^ o is the Larmor radius of the 
particles that dominate the cosmic ray number density at the shock, namely, for typical spectra 
of astrophysical interest, the lowest energy cosmic rays at the shock. At the shock location, the 
minimum momentum is the injection momentum, while at larger distances from the shock, the 
minimum momentum is determined by the diffusion properties upstream and is higher than the in- 
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jection momentum, since higher energy particles diffuse farther upstream. When the non-resonant 
mode exists, its maximum growth is found at krio » 1 . There may potentially be many impli- 
cations of this difference: the particle-wave interactions which are responsible for magnetic field 
amplification also result in particle scattering (diffusion). The diffusion properties for resonant and 
non-resonant interaction s are in general diffe r ent. Th e case of resonant interactions has been stud- 
ied in the literature (e.g. lLagage & CesarskyI (|1983al) ). at least for the situation 6B/B «; 1, but the 
diffusion coefficient for non-r esonant interactions (in either the linear or non-linear case) has not 



been calculated (see however 



Zirakashvili & PtuskinI (|2008|) ). The difference in wavelengths be- 



tween the two modes, in addition to different scattering properties, also suggests that the damping 
will occur through different mecha nisms. 

The calculation of IbcII ( 2004 ) has however raised some concerns due to the following three 
aspects: 1) the background plasma was treated in the MHD approximation; 2) a specific choice 
was made for the current established in the upstream plasma to compensate for the cosmic ray 
(positive) current; 3) the calculation was carried out in a reference frame at rest with the upstream 
plasma, where stationarity is in general not realized (although for small scale perturbations, the 
approximation of stationarity may be sometimes justified). 

In the present paper we derive the dispersion relation of the waves in a purely kinetic approach 
and investigate different scenarios for the microphysics that determines the compensating current. 
We show that the fastly growing non resonant mode appears when particle acceleration is very 
efficient, but whether it dominates over the well known resonant interaction between particles 
and alfven modes depends on the parameters that characterize the shock front, its Mach number 
pri marily . 



Belli (|2004|) also investigated the developement of the non-res onant modes by using nu merical 



MHD simulation s. His results have be en recently confirmed by 



Zirakashvili et al. 



similar approach. 



Niemiec et al. 



(120081) with a 



(|2008|) made a first attempt to investigate the development of the 



non-resonant modes by using PIC simulations. In this latter c ase, the auth ors find that the non 



resonant mode saturates at a much lower level than found by iBelll (|2004l) . However, as briefly 
discussed in § |5l these simulations use a set up that makes them difficult to compare directly with 
Bell's results. 

The paper is organized as follows: in § |2] we derive the dispersion relation of the unstable 
modes within a kinetic approach and adopting two different scenarios for the compensation of the 
cosmic ray current, namely compensation due to the motion of cold electrons alone (§ 12.11) . and 
to the relative drift of protons and electrons (§ 12.21) ; in § [3] we discuss the relative importance of 
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the resonant and non-resonant modes depending on the physical parameters of the system; we also 
derive analytic approximations for the large (§ 13.11) and small (§ 13.21) wavenumber limits; finally 
in § IHwe study the different modes during the Sedov evolution of a "typical" supernova remnant 
and for different assumptions on the background magnetic field strength; we conclude in § [51 

Throughout the paper we use the expressions accelerated particles and cosmic rays as referring 
to the same concept. 



2 THE KINETIC CALCULATION 

In this section we describe our kinetic calculation of the linear growth of waves excited by stream- 
ing cosmic rays upstream of a shock. This type of analysis is not suited for the computation of 
the saturation level of the instability, but it can be used to investigate the type of modes that may 
possibly grow to non linear levels. 

In the following, all calculations refer to the shock location and not to an arbitrary location 
in the plasma upstream of the shock. This is to say that the minimum momentum of the particles 
considered in our calculations is the injection momentum. 

In the reference frame of the upstream plasma the gas of cosmic rays moving with the shock 
appears as an ensemble of particles streaming at super- Alfvenic speed. This situation is expected to 
lead to streaming instability, as was indeed demonstrated in several previous papers (see 



Krall & Trivelpiece 



19731) for a technical discussion). 

In the reference frame of the shock, cosmic rays are approximately stationary and roughly 
isotropic. The upstream background plasma moves with a velocity towards the shock and is 
made of protons and electrons. The charge of cosmic rays, assumed to be all protons (positive 
charges) is compensated by processes which depend on the microphysics and need to be investi- 
gated accurately. 

The X-axis, perpendicular to the shock surface has been chosen to go from upstream infin- 
ity {x = -oo) to downstream infinity {x = -i-oo). Therefore a cosine of the pitch angle n = +\ 
corresponds to particles moving from upstream towards the shock. 

The dispersion relation of waves in thi s composite plasma, in the test-particle regime that we 



wish to investigate here, can be written as (|Krall & Trivelpiecd (|1973h ): 
c'^' 1 + ^""^ll r r' pMp)a -yu') \dfa , (kv . \ 1 df^ 



dp \u) p d/j. 



(1) 



oji^ ^ CO Jo J-1 CO + kv{p)ix + Q.a 

where the index a runs over the particle species in the plasma, a» is the wave frequency corre- 
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spending to the wavenumber k and Q.a, is the relativistic gyrofrequency of the particles of type a, 
which in terms of the particle cyclotron frequency f2* and Lorentz factor y is Qq, = Q* /y. For the 
background plasma and for any population of cold electrons one has Qa ~ ^*a- 

The positive electric charge of the accelerated cosmic rays, assumed here to be all protons, 
with total number density Ncr, must be compensated by a suitable number of electrons in the 
upstream plasma. In the following subsections we discuss two different ways of compensating 
the cosmic ray current and charge. In the first calculation we assume that there is a population 
of cold electrons which is at rest in the shock frame and drifts together with cosmic rays. These 



electrons cance 



Zweibell (11979. 



exac tly the positive charge of cosmic rays. This approach is similar to th at of 



Bell 



20031) and resembles more closely the assumptions of the MHD approach of 
(|2004|) . In the second calculation we assume that the current of cosmic ray protons is compensated 
by background electr ons and protons flowing at different speeds. This approach is similar to that 



of 



Achterberg 



mm- 



2.1 Model A: cold electrons 

Let n, and n^ be the number density of ions (protons) and electrons in the background plasma 
upstream of the shock. In this section we consider the case in which a population of cold electrons 
with density ncoid streams together with cosmic rays and compensates their charge. Therefore 
n^ = ni and n^oM = Ncr- In terms of distribution functions, the four components can be described 
as follows: 

fi(p,IJ^) = :^pSip-miV,)6(p-l) (2) 

feip,n) = :^^6ip -meV,)6(M- I) (3) 
2np^ 

fcR(p) = ^g(P)- (5) 

In the latter equation, which describes cosmic rays, g{p) is a function normalized so that f""" dp p^gip) 
1 . In the expressions above the background ions and electrons have been assumed to be cold (zero 
temperature). Introducing the thermal distribution of these particles does not add, as a first approx- 
imation, any important information to the analysis of the stability of the modes. One should check 
however that damping does not play an appreciable role, especially for the modes with high k (see 
the paper by Everett et al. (in preparation)). 
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The contribution of the background plasma of electrons and protons to the right hand side of 
Eq.[I]is easily calculated to be: 

Ane^Hi CO + kvs Ane^Ue co + kvs 



Cip-Mj CO + kVs ± O.* CO^nie CO + kVs ± 

Similarly the cold electrons with density Ncr contribute a term: 
Ane^ Nr 



(6) 



(7) 



CO irieico + Q*) 

The calculation of the cosmic ray contribution is slightly more complex. In its most general form, 
it can be written as 

ne^NcR r , idg , ^ ' 



XCR = 



dpv(p)p — d/u 
Jo dp J_i 



CO Jo ap j_i CO + kv{p)fi ± Qj 

The integral in the variable /i is now 



dfi ^, , , = r dn ^ - in d/u (1 -//) 6(kvn ± U 



), 



(8) 



(9) 



where P denotes the principal part of the integral and we have neglected co with respect to 
(low frequency modes). Using Plemelj's formula for the first term one obtains for the cosmic ray 
response: 

ne^NcR r~ , dg\^ 2 .,x2\, ^+p/Pm 



XCR 



cok 



1 + p/p„ 



^PminP 



n^e^NcR 



-i- 



cok 



dp -y- [p^ - Pminikf) , 
P,„M dp ^ ' 



where we have introduced the minimum momentum Pminik) = miQ.*/k, which comes from the 
condition that the second integral in Eq.[8]is non vanishing only when \fi\ < 1, namely when 



Qj Q.* Q.* 

V(p) ^ -r = ^ P = yf^iVip) ^f^i^= Pmm(k). 

K K'Y 



(10) 



The physical meaning of pmin is that of minimum momentum of the protons that can have a reso- 
nant interaction with waves of given wavelength. 

In the limit of low frequencies that we are interested in, co + kv^ O.* «; the contri- 
bution of the background plasma can be Taylor expanded and the unity in the dispersion relation 
(displacement current) neglected. So the dispersion relation reads 
Ncr. 



vlk^ = cb^ + 



-{cb - kvM [1 ± If{k) + ih{k)] , 



(11) 



where va = Bq/ ^jAmnini is the Alfven speed, co = co + kv^ is the wave frequency in the reference 
frame of the upstream plasma and we have introduced 



h(k) = — ^ — J dp— [p -Pminik) jln 



1 + p/p„ 



1 + P/P„ 



^PminP 



(12) 
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h(k) 



(k) 



dp—(p^-Pmm(kf)- 

P„u.m dp ^ ' 



(13) 



One should notice that the phase velocity of the waves in the plasma frame is = co/k and we 
want to concentrate on waves which have a velocity much smaller than the fluid velocity v.v (which 
is supersonic), therefore io kv^. In this limit, and using the fact that If = +If, one can write the 
dispersion relation as 

NcR-, 



'-kv,£i*{\^l\(k) + ih{k)\. 



(14) 



Befl 



(|2004|) by treating the background 



This dispersion equation is the same as that found by 
plasma in a MHD approximation. Here we have obtained Eq.[T4]by assuming that the cosmic ray 
current and charge are compensated by a population of cold electrons r noving with the cosmic rays. 



Bell (2004. At least within 



a setup which is equivalent to that of having the cosmic ray current of [ 
the context of this specific choice of the c ompensating current, treating the background plasma 
within the MHD approximation as done by iBelll (|2004|) does not change the results. However the 
question arises of whether the resulting dispersion relation may be different for a different and 
equally reasonable choice of the compensating current. In order to investigate this issue, in the 
section below we study the case in which the cosmic ray current is compensated by assuming a 
slow drift between thermal i ons and elec t rons. This approach resembles more closely the kinetic 



approach first put forward by lAchterbergl (Il983h. It is also importa nt to stress that this recipe is the 



same recently used in the PIC simulations of 



Niemiec et al. 



( 2008 ). 



2.2 Model B: compensation by electron-proton relative drift motion 



The approach described in this section is the one originally put forward by lAchterbergl (| 1 9 8 31) . We 
show that the dispersion relation is identical to that found in the previous section, provided that 
the density of cosmic rays is low enough compared with the density of the gas in the background 
plasma. 

Within this approach the electric charge of cosmic rays (assumed to be all protons) is compen- 
sated by the charges of electrons and protons in the background plasma 

A^CR + ni = ne, (15) 

and the total current induced in the background plasma by the presence of cosmic rays vanishes, 
namely 

= «,Vv - neVe- (16) 
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This condition can be realized by requiring that electrons and protons move with slightly different 
velocities, and respectively. The small drift between the two species leads to a current which 
compensates the cosmic ray current. 

For a given cosmic ray number density Ncr the contribution of accelerated particles to the 
dispersion relation does not change compared with the previous model of current compensation. 
The main differences with respect to the case presented in the previous section are that there are no 
cold electrons, and that electrons and protons in the background plasma have different velocities 
and densities. The contribution of the background plasma to the dispersion relation is then: 



Ane^rii oj + kvs Ane^rie co + kvg ^^^^^ 
(jcP-nii u) + kVs ± f2* (xP-nie o) + kvg + 

Now, introducing the frequencies o), = ct> + kVs and oje = co + kve, and using Eqs.[T5]and[T6l and 
again taking the low frequency limit, we obtain 

Ane^ I HjC I a)i \^ ni _ _ n,c _ NcrC 



^ ,+a)i — + Cbe + O^e >. (18) 



Now we notice that 



+ a>e = ±kvs— ~ ±kVs- 



Ncr + Ui Hi 

so that, after neglecting terms 0((NcR/ni)^y the contribution of the background plasma to the 
dispersion relation becomes: 

±]'^^Jl]\^]\ (19) 

VaI m OJ \VaI \oj/ 

It is easy to recognize that the first term is the same as the contribution of the background plasma 
in Model A, while the second term is equal to the contribution of the cold electrons in Model A. It 
follows that the two dispersion relations are identical up to linear terms in series of (NcR/ni) and 
start to differ at order O {{NcRlnif^. 



At this point it is worth pointing out that in the numerical PIC simulations of 



Niemiec et al. 



(|2008|) the compensating current is realized by assuming a drift between protons and electrons, 
exactly as discussed in this section. However, in order to be able to carry out the calculations, the 
authors are forced to adopt unrealistically large values of the ratio NcR/rii (for the most realistic 
cases they use Ncr/iii = 0.3), which, as discussed above, do not necessa rily result in a dispersion 



relation with the same characteristics as that found above and derived by 



Belli ( 20041) 
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3 RESONANT AND NON RESONANT MODES 

In this section we investigate the modes that result from the dispersion relation in Eq. [141 For the 
sake of simplicity we carry out our calculations for a power law spectrum of accelerated particles 
with the canonical shape g(p) oc p""*, which is expected from diffusive acceleration at strong 
shocks in the test particle regime. This regime may potentially be different from the one relevant 
for particle acceleration in supernova remnants, in which acceleration takes place efficiently and 
the shock becomes cosmic ray modified by non-linear effects. The most evident effect of this 
non-linear reaction is the formation of a shock precursor which reflects in concave spectra of 
accelerated particles. The energy content of these accelerated particles is typically dominated by 
the highest energy particles while numerically the low energy particles are still dominant. On the 
other hand, efficie nt acceleration must p roceed parallel to effective magnetic field amplification, 
and as showed by ICaprioli et al.l (l2008allb() . the dynamical reaction of the strong field leads to a 
smoothening of the shock precursor, which in turn leads to spectra of accelerated particles that 
though still concave, are closer to power laws. It is straightforward to grasp the deeply non-linear 
nature of the problem at hand, where one is trying to describe the growth of waves associated to 
cosmic ray streaming, but the latter only occurs after the waves have grown to interesting levels, 
so as to guarantee efficient particle acceleration. At the present time a full analysis of the whole 
problem is simply not achievable. In the following we still adopt the assumption of power law 
spectra of accelerated particles, and we point out, where necessary, if a result may be substantially 
modified by non-linear effects. 

The suitably normalized distribution function g{p) is 



= 11 P. 



-4 



g(p) = ^ — ®ip-Po) ®{Pmax - P) (20) 

pI \PoI 

where is the step function and takes into account the limited range of momenta spanned by the 
cosmic ray particles in the acceleration region. 

The momentum p^ is the minimum momentum of accelerated particles at the shock location, 
namely the injection momentum. Below, we discuss in detail the dependence of our results on the 
choice of the value of po- 

Let us now consider the integrals h and h in Eqs. [T2lfT3l We integrate by parts after performing 
the substitution 

- fcr^,o, (21) 



Pminik) PO 

obtaining: 
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1 



nik) = -r 



Pi 



h(k) 





\ + s 


+ 2*1 


g{s){{s' -l)\n 


1-5 





Jo 



ds s g{s) In 



1 + s 



1 



dssg(s) 



(22) 
(23) 



After defining sq = kr^o and S2 = krifiPmax/Po and using the expression for ^(5) obtained from 
Eq.l20lone finds: 



2I0 



ds s In 



l + s 



So 

'4 



i|(.^-l)ln 



l + s 



\-s 



-2 s 



so 



n 



h(k) = 2 



with si = Max[l, sq]- Finally: 



i'2 



It(k) 



1 



4-krL,o 



(ikrtflf - 1) In 



1 + krLfl 



- 2 fcr^j 



n 



1 - krtfl 
krtfl krifi < 1 
ikrL,or' krL,o > 1 ■ 
In terms of the latter, the imaginary and real parts of the frequency can be written as: 



h(k) = -- 



^]ik) = ]- 



COR{k) = - 



- (k^vl + a{l+ h(k))) + ^J{k^vl ±a{l+ h(k))) + a^L 

ah 
2u)i 



(24) 
(25) 

(26) 
(27) 

(28) 
(29) 



where a = -^fcViil*. It is useful to express a as a function of the acceleration efficiency of the 
shock. The total pressure in the form of accelerated particles is 

Pmax 



-NcR dpp^v(p)g{p) ~ -NcRcpo In 
J PI 



'pa 



niiC 



(30) 



The second part of this expression is a consequence of the fact that for spectra harder than p^^ 
and not harder than p''^, which are of interest here, the pressure is mainly contributed by relativis- 
tic particles with p ~ m,c. In case of cosmic ray modified shocks, spectra can becomes mildly 
harder than p""*, but this is not expected to aff"ect our conclusion in a dramatic way, since it only 
introduces a weak dependence on Pmax- Moreover, in case of strongly modified shocks, the major 
complication does not come from the spectrum of accelerated particles, but rather from the fact 
that the upstream plasma develops a precursor, namely a gradient in the fluid velocity that makes 
the standard treatment illustrated here formally not applicable. On the other hand, as pointed out 
several times throughout the paper, the magnetic reaction of the shock acts in the direction of 
smoothening the precursor, so that although formally the mathematical treatment illustrated here 
is not applicable, in practice it should provide a good description of the physical processes at work. 



> 



> 



A kinetic approach to growth rates of streaming instability at supernova shocks 1 1 



10^ 
10 

1 

10"^ y- 



10" 
10" 



10" 



non-resonant mode 




resonant mode 



Re(cj) 



Im(cj} 



10"^ 



10 



10^ 



10^ 



10^ 



10' 



krL,o 



Figure 1. We plot the real and imaginary part of the frequency as a function of wavenumber for the resonant and non-resonant modes. Wavenumbers 
are in units of l/r^o, while frequencies are in units of Kcrt o). The top panel refers to the non-resonant branch, while the lower panel is for the 
resonant branch. In each panel, the solid (dashed) curve represents the real (imaginary) part of the frequency. The values of the parameters are as 
follows: V5 = 10' cm s"'. Bo = l/JC?, «, = 1 cm"'', = 0.\,p„ax = 10^ mpC. 



If we define 77 = Pdininiivl) as the acceleration efficiency, we can write: 

^ Iv] k k 

a = 3t] — ■ = cr — , 

R c r,Q rro 



(31) 



where R = In(^j^) and r^ o = poc/eBo is the Larmor radius of the particles with momentum po 



in the background magnetic field Bq. We have also introduced cr = 3r]-^ — . A resonant mode can 



be obtained from Eq. [28] with both signs of the polarization. On the other hand the non-resonant 
mode only appears when the lower sign is chosen. 
We notice that the following relation holds: 



_ NcrPq_VsC_ _ An ^rifl 
v\ ni nijC v\ c Bq' 



(32) 



where J = eNcuVs- This means that the system is strongly current driven when » 1. 

Therefore the parameter cr/v^ controls the growth rate of the non-resonant mode: when cr/v^ » 
1 the non-resonant mode is almost purely growing and its growth is very fast. When cr/v^ «; 1, 
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the non-resonant mode is subdominant and a resonant mode is obtained, asymptotically identical 
to that corresponding to left-hand polarized waves (upper sign in Eq.[28l). 

In the following we often refer to the mode arising with the lower sign of the polarization as 
the non-resonant mode, although one should keep in mind that its peak growth rate reduces to that 
of the standard resonant mode in the limit cr/v^ «; 1. 

In Fig. [H we plot the solution of the dispersion relation in a case for which cr/v^ » 1. The 
values of the parameters are: Vs = lO^cm s \ Bq = l/uG, ni = 1 cm"^, 77 = 0.1, p^ax = 10^ i^pC- 
The frequency (y-axis) has been normalizated to the advection time for a fluid element upstream 
of the shock through the characteristic distance cri o/vs, namely cr^o/Vj . It is worth stressing 
however that this is a good estimate of the diffusion time scale only in the case of Bohm difl'usion, 
when the difl^usion coefficient is D{p) « nflC and the diffusion time scale is D{p)/vs « cr^ o/v^. 
The plots in the upper (lower) panel in Fig. [T] are obtained by choosing the lower (upper) sign of 
the polarization in the dispersion relation (Eq. [28]and Eq.[29l). 

First, let us comment on the consistency of our derivation of the disperion relation. It is easy to 
check, from Fig.[Tl that these are indeed low frequency modes. More specifically, they satisfy both 
assumptions underlying our calculation: co «c kvs and o) <^ Qj. Moreover, the non-resonant mode 
(lower sign of the polarization in the dispersion relation) is characterized by an imaginary part that 
is much larger than its oscillatory part for a very large range of wavenumbers. In this same range 
of k, for our choice of the parameters, its growth is much faster than for the resonant branch. 

Further insight in the behaviour of the different wave modes can be gained by investigating 
the limits of the dispersion relation for the regimes krio 1 and kr^ o » 1. Based on these 
asymptotic trends it is easy to explain (see ^3.31) why the normalized frequencies as functions of 
the normalized wavenumber, krifl, do not depend on the choice of the minimum momentum. 



3.1 Large wavenumber limit: kr^) » 1 

For krio » 1 one easily obtains that Ii(k) =^ -(l/3)/(fcri_o)^ and l2{k) - -(n/4)/krLfi. In Eq.l28] 
there are three terms that determine the actual dependence of a>/ on wavenumber k, namely 

k OC k , a(l + I\) ~ cr — oc k, ah ~ -7-5— °^ k . 

In the range kr^fi » 1 the third term is always subdominant. Moreover we can identify two critical 
values of the wavenumber k, ki and k2, such that for » » 1 / the first term dominates upon 
the third, and for ^ » ^2 the first term also dominates upon the second term, which is linear in k. 
It is easy to find that 
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\lyuG/ Icm- 



Similarly 

For k <^ k2 Eq. 128) gives 

a;,^ = icrALi + l + ^ 1 1 (35) 

When the upper sign is chosen in the above equation, we obtain wj oc k'^^^, while when the lower 
sign of the polarization is considered one finds cbj oc fc'^ ^, namely t he growth rate of the waves 



Belli (120041) . For this mode cbj increases 



increases with k. This is the non resonant branch found by 
with ^ up to ~ ^2 and in the range of wavenumbers between ki and ^2 is larger than for the 
resonant waves. The maximum growth rate is obtained for k ~ ^2- In fact, for k kx one finds: 

64 v^r^Q (krioY 

which implies d)j oc k~^ for both the resonant and non-resonant modes. 

The non resonant mode disappears when ki becomes larger than kj, which happens for 

i < i - < « X -i'"(T^))(lo^[' (l^f (^)" • 

For the reference values of the parameters, the non-resonant mode grows faster than the res- 
onant mode only for unreasonably low efficiencies of particle acceleration, as one may conclude 
by comparing Eq. [37] with ?; ~ 0.1 - 0.2 required for the association of cosmic rays to supernova 
remnants. 

On the other hand, for shock velocity = lO^cm s ' and magnetic field Bq = 5//G one easily 
sees that the limit in Eq. [37] becomes rj < 0.1. This implies that the resonant and non resonant 
modes compete during the history of a supernova remnant, with the resonant mode prevailing 
during the stages in which the shock has slowed down appreciably. We will comment further on 
this point below. 

We summarize the results of this section by giving the following useful approximations to the 
solution of the dispersion relation in the large wavenumber limit. For l/ri o «c «; ^ «; ^2 we 
have: 




ojY' « aj~'' ~ J k"^''- and o^f""''' « ojT « -./^ k^'^ (38) 
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For A:2 «; we have: 



-res 



~, non-res 
0)1 



n cr 



and 



VAr 



-res 



~, non-res 



kVA 



(39) 



L,0 



3.2 Small wavenumber limit: kri Q «; 1 

In the limit of perturbations with wavelength much larger than the gyroradius of the lowest energy 
particles in the cosmic ray spectrum, the results depend again on the ratio between v\ and cr. As 
we already mentioned, for most regions of the parameters space v\ «; cr, to which case we limit 
our analysis here, while we defer to next section a discussion of what happens for slow shocks. 



For kriQ 1 one has Ii{k) -1. Hence, in the limit v\ «; cr, we find for both signs of 
polarization: 




3.3 Dependence on the minimum momentum 

The number density of accelerated particles at the shock depends on the minimum momentum 
Pq, which in practice it is hard to know or to predict, since it is determined by details of the 
microphysics of the shock formation. It is therefore useful to address the issue of the dependence 
of our results from the assumed value of po- 

Let us start from the non resonant mode. For kr^Q » 1, following Eq. [35l we have that 
ib] oc kjrio for k k2 and oc 1 /(k^rj^^) for k :s> k2.It follows that apparently the dependence 
on po is severe. However, if the wavenumbers k are normalized to l/ri o and <&/ is normalized to 
the diffusion time scale rtoc/vl, the resulting imaginary and real parts of the frequency are inde- 
pendent of the choice of pq. In other words, the absolute value of the growth rate is affected by the 
choice of po, but the rate itself, which is regulated by the time available for the instability to grow, 
as measured by the time scale riflc/vj, is independent of the choice of pq. Similar considerations 
can be easily repeated for the resonant mode and for the low frequency regime of the non resonant 
mode. 



4 RESONANT AND NON-RESONANT MODES IN SNRS 

We now study the relative importance of the resonant and non-resonant wave modes during the 
evolution of a SNR. We consider a remnant originating in a SN explosion with energy Esn- Once 
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the remnant has entered the Sedov phase, the shock velocity as a function of time t can be written 
as: 

\l05ierg/ llcm-V \lOV/ 
As discussed in the previous section, the existence of the non-resonant mode depends on the ratio 
cr/v^, which can be written as a function of the age of the remnant as: 

This implies that during the evolution of a "typical" supernova remnant. Bell's instability, 
which requires cr/v\ > 1, is likely to operate only at early times after the beginning of the Sedov 
phase. The non-resonant mode disappears when the remnant is a few 10"^ yr old if it is expanding in 
a ljuG magnetic field and 10 times faster (age about a few lO"* yr) if the background magnetic field 
is 10 times higher. At later times, the streaming cosmic rays will still amplify the field but only via 
the classical resonant mechanism. This is also clear from Fig.[2]where we plot the growth rate of 
the non-resonant mode as a function of age for the above mentioned values of the magnetic field: 
Bo = 1/iG in the upper panel and Bq = lOfxG in the lower panel. From the plots in Fig. El where 
again the time-scale for wave growth is normalized to the fastest of the time-scales involved in the 
system dynamics, cr^ o/v^, one immediately sees that at least the resonant mode of the streaming 
instability still grows efficiently after 10^ yr since the supernova explosion. The non-resonant 
mode, on the other hand, soon becomes subdominant. 

The non-resonant mode grows the fastest at ^ ~ ^2/2, so that from Eq. [38] we can derive the 
maximum growth rate as: 

.1/2/7 \l/2 

cr \ kj\ cr 



Tmax ~ maxioji) = — — = — . (43) 

{rtfll \2/ 2^'^VArL.o 

It is clear that, remarkably, r^ax does not depend on the background magnetic field Bq. On the other 
hand, the wavenumber at which the growth is maximum does depend on Bq (see Eq. U^ . These 
trends are clearly seen from Fig. [31 where we plot the dependence on time of the maximum growth 
rate, Fmax, and of the wavenumber for which this occurs, for both the resonant and non-resonant 
modes. The plot refers to the "typical" SNR parameters considered above and the two mentioned 
values of the background magnetic field strength. In both this figure and Fig. [2] a 10 % particle 
acceleration efficiency was assumed, and kept constant during the evolution of the remnant. This 
latter assumption is definitely not very realistic, as can be demonstrated by using non-linear theory 
of particle acceleration. 
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Figure 2. We plot the growth rate of the non-resonant mode as a function of wavenumber. Wavenumbers are in units of l/r^ o, while growth-rates 
are in units of the advection time v^, /cr^ o- The different curves in each panel refer to different ages of the remnant: solid is for 10^^ yr, dashed is 
for 5 X lO', dot-dashed is for lO*, dot-dot-dashed is for 5 X 10^ and finally dotted is for 10^. Two dift'erent values of the background magnetic field 
strangth are assumed in the two panels: Bo = 1 I^G in the top panel, and Bq = 10/jG in the bottom one. The shock velocity is computed according 
to the Sedov expansion of a remnant with Eg^ = 10^' erg. The remaining parameters are as follows: n, = 1 cm"^, rj = 0.\, Pi„ax = lO' nipC. 



In Fig. |3]the time at which the fastest growing mode switches from non-resonant to resonant is 
identified by the intersection between the dashed line and the solid (Bq = 10/iG) or the dot-dashed 
(Bq = 1/iG) one depending on the magnetic field strength. The dominant wave mode progressively 
moves to larger wavelengths. The implications of this peculiar trend are expected to be profound 
on the determination of the diffusion coefficient: we recall that the standard Bohm diffusion is the 
limit obtained for resonant interactions of particles and waves when 6B(k) = Bq for any value of 
k. For non-resonant modes, the diffusion properties need to be recalculated from first principles. 
On one hand, since the most unstable modes have k ^ \ /rio, most particles do not resonate 
with these modes and the typical deflection suffered by a single particle within a spatial scale 
~ l/k is very small. On the other hand the number of scattering events is verv large, therefo re a 



substantial reduction of the di ffusion coefficient can still be expected (see 
iZirakashvifi & Ptuskinlfcoosl) ). 



Reville et al. 



E008)and 
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Figure 3. In the top panel we plot the maximum growth rate of the resonant and non-resonant branches as a function of the age of the supernova 
remnant. The growth rate Fmax = max(i3j;) is in units of yr"' while along the x-coordinate time is expressed in yr. The notation for the different 
curves is as follows: the dashed line and the symbols are for the non-resonant mode in a 10/iG and IpG magnetic field respectively; the solid 
and dot-dashed lines are for the growth of the resonant mode, again for Bq = 10/iG and Bg = l^iG respectively. In the bottom panel we plot the 
wavenumber corresponding to the fastest growing wave mode for the same situations considered above. Wavenumbers are in units of r^.o and the 
notation for the different line-types is as follows: the dashed and dotted lines are for the non-resonant mode in a lOyuG and IfiG magnetic field 
respectively; the solid line and symbols are for the resonant mode, again for Bo = 10/iG and Bq = ^^^G respectively. The shock velocity changes 
with time according to the Sedov evolution of a remnant with Egn = 10^' erg. The remaining parameters are as follows: n, = 1 cm"', = 0.1, 

Pmax = 10^ m,,C, pi„j = lO^^m/C. 



5 CONCLUSIONS 



We investigated the excitation of streaming instability induced by accelerated particles in the 
vicinity of a non-relativistic shock wave, typical of supernova shells expanding in the interstellar 
medium. The calculation is based on kinetic theory, hence we do not require the MHD approxima- 
tion to hold for the background plasma. We find that the dispersion relation of the waves leads to 
the appearance of two mo des, a resonant a nd a non-resonant one. T he former is the well known un- 



stable mode, discussed by 



Zweibell ( 1979 ) and 



AchterbergI (|1983|) . bas ed on a reson ant interaction 



between waves and particles. The latter is similar to that discussed by 



Belli (12004 . who however 



based his a nalysis on a set of assumptions that called for further investigation: the calculation of 



Bell 



2004() is based on the assumption that the background plasma can be treated in the MHD 



regime, and makes specific prescriptions on the return current which compensates the cosmic ray 
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current upstream of the shock. Moreover, the whole calculation is carried out in the frame of the 
upstream plasma, where in principle there is no stationary solution of the problem. 

Our kinetic calculations are carried out for two models of the compensating current: in the first 
model, the return current is established through a population of cold electrons, at rest in the shock 
frame, which exactly compensate the positive charge of cosmic ray protons. In the second model, 
the return current is due to a slight drift between ions and electrons in the background plasma 
upstream of the shock. We have demonstrated that the dispersion relation of the waves is the same 
in the two cases, to order C)(A^cr/"i)^- 

The resonant and the non-resonant mode are found at the same time, with growth rates which 
in the general case are different. The non-resonant mode is almost purely growing and is very 
apparent when particle acceleration is efficient. The parameter that regulates the appearance of 
the non-resonant mode is cr/v^, where cr = 3t]vI/(cR). When cr/v^ » 1, the waves excited in 
a non-resonant way grow faster than the resonant modes and may lead to a substantial magnetic 
field amplification. 

The strong dependence of cr on the shock velocity implies that the non-resonant mode is likely 
to be the dominant channel of magnetic field amplification in SNRs in the free expansion phase 
and at the early stages of the Sedov-Taylor phase of adiabatic expansion. At later times, the non- 
resonant mode collapses on the resonant mode, which keeps providing appreciable growth for 
longer times, at least if damping mechanisms are neglected. The growth of the fastest non-resonant 
mode is independent on the strength of the unperturbed initial magnetic field Bq. 

The non-resonant mode, when present, grows the fastest at wavenumber ^2 given by Eq. [341 
which in the cases of interest is much larger than l/rio, where r^ o is the gyroradius of the par- 
ticles with minimum momentum in the cosmic ray spectrum. These modes are therefore short 
wavelength waves, which is the m ain reason why the assumption of stationarity in the upstream 



frame, as required by 



Bell 



(|2004f) . was acceptable, despite the impossibility of reaching actual 
stationarity in that frame. 

The numerical results in this paper were specialized to the case of a power law spectrum 
p~'^ of accelerated particles, typical of Fermi acceleration at strong shocks. However, one should 
keep in mind that the levels of efficiency required for the non-resonant mode to appear are such 
that the dynamical re a ction of the accelerated particles on the shock cannot be neglected (see 



Malkov & Q'C Druryl (|200lh for a review). This backreaction leads to several important effects: 



on one hand the spectra of accelerated particles become concave, and concentrate the bulk of 
the energy in the form of accelerated particles at the maximum momentum. On the other hand, the 
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efficiently amplified magnetic field also exerts a strong dynamical r eaction on the system, provided 



Caprioli et al. 



2008a bh). This 



the magnetic pressure exceeds the gas pressure in the shock region (l< 
second effect results in an enhanced acceleration efficiency (due to large B-fields) but weaker shock 
modification (spectra closer to power laws) due to the reduced compressibility of the plasma in the 
presence of the amplified magnetic field. These non-linear effects cannot be taken into account in 
the type of calculations presented here, although we do not expect the qualitative character of our 
conclusions to be affected profoundly by them. 

Nevertheless it is useful to go through the possible consequences of the non-linear effects in 
a somewhat deeper detail: one can expect two types of complications, one of principle and the 
other in the numerical values of the growth rates. The latter simply derives from the approxima- 
tions intrinsic in the assumptions we made: for instance the adoption of a power law spectrum of 
accelerated particles with slope 4, which clearly fails when a precursor is formed. 

There are then complications coming from deeper unknown pieces of Physics or from a not to- 
tally satisfactory mathematical approach. For instance, the standard perturbative approach adopted 
here is based on the assumption of a spatially uniform background. The presence of a precursor 
invalidates this assumption, although probably not in a dramatic way. 

Since the non-resonant mode appears for large values of k, the relevant quantities can be as- 
sumed to be spatially constant in the precursor on scales ~ 1 /k, so that in this respect our calcu- 
lations are still expected to hold, and probably to a better accuracy for the non-resonant modes 
{k » l/^L,o) than for the resonant ones {k ~ I/^l.o)- Moreover, as stressed above, the dynamical 
reaction of the magnetic field leads to weaker modification of the shock, and therefore to spec- 
tra with less prominent concavity (closer to p''^). Also in this respect, the calculations presented 
here should serve as a good description of all relevant physical effects related to the growth of the 
cosmic ray induced instabilities. 

More important, the acceleration process is directly affected by the physics of particles' dif- 
fusion in the shock region, which in turn is determined by the excited waves. This intrinsic non- 
linearity cannot be taken into account in perturbative approaches like ours or like Bell's, and one 
should always be aware of this limitation. Even more, while the diffusion coefficient for resonant 
modes can at least be derived in quasi-linear theory, at present there is n o derivation of the dif- 



Zirakashvili & Ptuskin 



fusion coefficient associated with scattering on non-resonant modes (see 
J2OO8 ) for a first attempt at discussing this effect). 

Another issue that deserves further investigation is that of determining the level of field am- 
plification at which the instability saturates. This cannot be worked out within a linear theory 
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calculation and onl y numerical simu 



have been made by 



Niemiec et al 



Belli (120041) and 



ations can address this issue. Recent efforts in this direction 



Zirakashvili et al. 



(|2008h through MHD simulations and by 



(|2008h by using PIC simulations. While the first two papers find a saturation level 



6B^ liAn) ~ (y Jc)Pc, in the third paper a much lower level of field amplification is found. The 
authors conclude that the existence of large magnetic field amplification through the excitation of 
non-resonant modes is yet to be established. 

Although we agree with this conclusion, we also think that the setup of the PIC simulation 



by 



Niemiec et al 



(|2008f ) is hardly applicable to investigate the excitation of the Bell instability at 
shocks, or at least several aspects of it should be studied more carefully. First, they carried out the 
calculations in a regime in which the condit ion of strong magnet ization, co «; ^^o^ was violated. 



Second, in order to carry out the calculations. 



Niemiec et al. 



(20081) are forced to assume unrealisti- 



cally large val ues for the ratio NcR/n j (of order 0.3 for their most realistic runs). The return current 



as assumed by 



same dispersion relation as 



Niemiec et al.l (I2008h co rresponds to our second model, which however leads to the 
Belli (|2004r ) only at order O {NcRlntf, which is not necessarily the case 



here. Moreover, the spectrum of accelerated particles is assumed to be a delta-function at Lorentz 
factor 2, instead of a power law (or more generally a broad) spectrum. It is not obvious that for 
non-resonant modes this assumption is reasonable. But the most serious limitation of this PIC sim- 
ulation is in the fact that the authors do not provide a continuous replenishment of the cosmic ray 
current, which is instead depleted because of the coupling with waves. In the authors' view this 
seems to be a positive aspect of their calculations, missed by other approaches, but in actuality 
the cosmic ray current is indeed expected to be stationary upstream, and we think that the PIC 
simulation would show this too if particles were allowed to be accelerated in the simulation box 
instead of being only advected and excite waves. Clearly if this were done, the spectrum of acceler- 
ated particles would not keep its delta-function shape, but should rather turn into a power-law-like 
spectrum. The latter issue adds to the absence of a replenishment of the current, which seems to 
us to be the main sho rtcoming of these simulations. Overall, it appears that the setup adopted by 



Niemiec et al 



(|2008|) would apply more easily to the propagation of cosmic rays rather than to 
particle acceleration in the vicinity of a shock front. 

The issue of efficient magnetic field amplification, possibly induced by cosmic rays, has be- 
come a subject of very active debate after the recent evidence of large magnetic fields in several 
shell-type SNRs. The implications of such fields for particle acceleration to the knee region, as 
well as for the explanation of the multifrequency observations of SNRs, are being investigated. 
Probably the main source of uncertainty in addressing these issues is the role of damping of the 
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excited waves. For resonant modes, ion-neutral damping and non-linear Landau damping have 
been studied in some detail: their role depends on the temperature of the upstream plasma and 
on the shock velocity. For non-resonant modes, being at high k, other damping channels could 
be important (Everett et al., in preparation). Whether SNRs can be the source of galactic cosmic 
rays depends in a complex way on the interplay between magnetic field amplification, damping, 
particle scattering and acceleration, together with the evolution of the remnant itself. 

ACKNOWLEDGMENTS 

The authors are very grateful to J. Everett and E. Zweibel for a critical reading of the manuscript 
and for ongoing collaboration and to B. Reville for constructive discussions. They also thank J. 
Kirk for identifying a problem in a previous version of the manuscript in the low-k part of the 
dispersion relation. This work was partially supported by PRIN-MIUR 2006, by ASI through 
contract ASI-INAF 1/088/06/0 and (for PB) by the US DOE and by NASA grant NAG5-10842. 
Fermilab is operated by Fermi Research Alliance, LLC under Contract No. DE-AC02-07CH1 1359 
with the United States DOE. 

REFERENCES 

Achterberg A., 1983, A&A, 119, 274 
Bell A. R., 1978, MNRAS, 182, 147 
Bell A. R., 2004, MNRAS, 353, 550 

CaprioU D., Blasi R, Amato E., Vietri M., 2008a, ApJ Lett, 679, L139 

Caprioli D., Blasi R, Amato E., Vietri M., 2008b, "Dynamical feedback of self-amplified mag- 
netic fields in cosmic rays modified shocks". Submitted to MNRAS, preprint arXiv:0804.2884 

Krall N. A., Trivelpiece A. W., 1973, Principles of Plasma Physics. McGraw-Hill, pp 249-257 

Lagage P O., Cesarsky C. J., 1983a, A&A, 118, 223 

Lagage P O., Cesarsky C. J., 1983b, A&A, 125, 249 

Malkov M. A., O'C Drury L., 2001, Reports of Progress in Physics, 64, 429 

Morlino G., Amato E., Blasi P., 2008, "Gamma ray emission from SNR RX J1713.7-3946 and 
the origin of galactic cosmic rays". Submitted to MNRAS 

Niemiec J., Pohl M., Stroman T., 2008, ArXiv e-prints, 802 

Pohl M., Yan H., Lazarian A., 2005, ApJ Lett, 626, LlOl 

Reville B., O' Sullivan S., Duff"y P, Kirk J. G., 2008, MNRAS, 386, 509 



22 E. Amato and P. Blasi 
Skilling J., 1975, MNRAS, 173, 255 

Volk H. J., Berezhko E. G., Ksenofontov L. T., 2005, A&A, 433, 229 
Zirakashvili V. N., Ptuskin V. S., 2008, ArXiv e-prints, 801 
Zirakashvili V. N., Ptuskin V. S., Voelk H. J., 2008, ArXiv e-prints, 801 

Zweibel E. G., 1979, in Arons J., McKee C, Max C, eds. Particle Acceleration Mechanisms 
in Astrophysics Vol. 56 of American Institute of Physics Conference Series, Energetic particle 
trapping by Alfven wave instabilities, pp 319-328 

Zweibel E. G., 2003, ApJ, 587, 625 



